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The role of the top quark in completing the Standard Model quark sector is re- 
viewed, along with a discussion of production, decay, and theoretical restrictions 
on the top quark properties. Particular attention is paid to the top quark as a lab- 
oratory for perturbative QCD. As examples of the relevance of QCD corrections in 
the top quark sector, the calculation of e+e - — > it at next- to-leading-order QCD 
using the phase space slicing algorithm and the implications of a precision mea- 
surement of the top quark mass are discussed in detail. The associated production 
of a it pair and a Higgs boson in either e+e~ or hadronic collisions is presented at 
next-to-leading-order QCD and its importance for a measurement of the top quark 
Yukawa coupling emphasized. Implications of the heavy top quark mass for model 
builders are briefly examined, with the minimal supersymmetric Standard Model 
and topcolor discussed as specific examples. 



1. Introduction 

Long before its discovery in 1995, 1 the top quark was regarded as an essen- 
tial ingredient of the Standard Model (SM) of particle physics. Its existence, 
and many of its properties, are determined by requiring the consistency of 
the Standard Model. These requirements (discussed in Section 2) specify 
the couplings of the top quark to the SU (3) x SU(2) l x U(1)y gauge bosons 
and many of the top quark properties: 
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• Weak isospin partner of b quark: Tg = i 

• Color triplet 



spin- i 
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The top quark mass, which is measured at CDF and DO to be M t = 
174.3 ± 5.1 GeV, 2 is not predicted in the Standard Model, but is restricted 
by precision electroweak measurements. 3 ' 4 In order to confirm that the ob- 
served top quark is that predicted by the Standard Model, all of the prop- 
erties listed above must be experimentally verified by direct observation. 
Section 3 contains a discussion of the measurements of top quark properties 
at the Tevatron and surveys the improvements expected at future colliders. 
Recent reviews of top quark physics can be found in Ref. 5. 

In Section 4, we discuss the top quark as a laboratory for perturbative 
QCD. At the top quark mass scale, the strong coupling constant is small, 
a s {M t ) ~ 0.1, and so QCD effects involving the top quark are well behaved 
and we expect a perturbation series in a s to converge rapidly. The process 
e + e~ — ► ti provides an example of QCD effects in top quark production and 
the next-to-leading-ordcr (NLO) QCD corrections are described in detail 
using the phase space slicing (PSS) algorithm. A knowledge of these higher 
order QCD corrections is vital for extracting a precise value of the top 
quark mass from the threshold behavior of the e+e~ — > it cross section. As 
a further example of the role of QCD effects in the top quark sector, we 
consider the associated production of tih and the implications for measuring 
the top quark Yukawa coupling. 

Finally, in Section 5, we discuss the importance of the top quark mass 
in model building. Since the top quark is heavy, it is expected to play 
a special role in elucidating the source of fermion masses. We begin by 
discussing the significance of the top quark in supersymmetric models, and 
finish with a brief discussion of the top quark and models with dynamical 
symmetry breaking. 

2. Who needs a top quark? 

In this section we discuss some of the reasons why the top quark was be- 
lieved to exist even before its experimental discovery. These considerations 
fall into three general categories: theoretical consistency of the Standard 
Model gauge theory (anomaly cancellation), consistency of b quark mea- 
surements with SM predictions, and consistency of precision measurements 
with the SM. We then turn to a discussion of top quark production and 
decay mechanisms at the Tevatron and the LHC. 

The particles of the first generation of fermions, along with the Higgs 
doublet, are shown in Table 1, with their gauge quantum numbers. The 
third generation is assumed to follow the same pattern, with the left-handed 
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top and bottom quarks forming an SU(2)l doublet with hyperchargc, Y = 
i. The right-handed top and bottom quarks are SU(2)l singlets. Our 
normalization is such that: 

Qem=T 3 +Y (1) 

withT 3 = ±i. 



Tablcl. Fcrmions in the first genera- 
tion of the Standard Model and their 
SU(3) X SU(2) L X U(1) Y quantum num- 
bers. 



Field 


SU(3) 


5(7(2) 


(uL,d L ) 


3 


2 


UR, 


3 


1 


d>R, 


3 


1 




1 


2 




1 


1 


H 


1 


2 



U(lh 



2.1. Anomaly Cancellation 

The requirement of gauge anomaly cancellation 6 ' 7 puts restrictions on the 
couplings of the fermions to vector and axial gauge bosons, denoted here 
by and . The fermions of the Standard Model have couplings to the 
gauge bosons of the general form: 

L ~ g^T^^A^ + g v ^T a ^V a ^ (2) 

where T a is the gauge generator in the adjunct representation. These 
fermion-gauge boson couplings contribute to triangle graphs of the form 
shown in Fig. 1. The triangle graphs diverge at high energy, 

/d n h 1 
(^p. (3) 

where r\i = =f1 for left- and right-handed fermions, tpL,R = 5(1 T 75)^- 
This divergence is independent of the fermion mass and depends only on 
the fermion couplings to the gauge bosons. Such divergences cannot exist 
in a physical theory, and must somehow be cancelled. The theory can 
be anomaly free in a vector-like model where the left- and right-handed 
particles have identical couplings to gauge bosons and the contribution to 



Figurel. Generic Feynman diagram contributing to gauge anomalies. The fermion loop 
contains all fermions transforming under the gauge symmetry, a, b, c are gauge indices. 



Eq. 3 cancels for each pair of particles. From Table 1, however, it is 
clear that the Standard Model is not vector-like. The anomaly, T abc , must 
therefore be cancelled by a judicious choice of fermion representations under 
the various gauge groups. 

The only non-vanishing contribution to the anomaly in the Standard 
Model is from 

-£Tr[Y{T a ,T b }}, (4) 

where T a are the SU(2)l generators and the sum is over all fermions in the 
theory. Eq. 4 vanishes for the hypercharge assignments given in Table 1. 
Note that the anomaly cancels separately for each generations of fermions. 
The cancellation of the gauge anomalies in the Standard Model for the third 
generation therefore requires that the b quark have a T 3 = i partner with 
electric charge Q\ m = | | e |, and hypercharge Y l — Q' m — T|. The partner 
of the b quark is by definition the top quark. Since anomaly cancellation is 
independent of mass, a priori , the top quark mass could be anything. 



2.2. b quark properties 



Many of the experimental properties of the b quark require that it be a 
-| particle with Q h em = — | and Y b = i. The coupling of the 
b quark to the Z boson can be tested to check if it has these quantum 
numbers. The SM fermions couple to the Z boson as, 
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where 

Ri = —2Qi sin 2 9w 

Li = 2T 3t - 2Qi sin 2 6 W , (6) 

and 9w is the electroweak mixing angle. The experimental value of Rhad 
is sensitive to Q b em , 

_ er(e + e~ — > hadrons) 

Rhad = , , — T _s ■ (') 

cr(e + e [V [i ) 

At a center-of-mass energy, -^s = 2m bl i?6 depends sensitively on the b 
quark electric charge, 

SR had (2m b ) = N c {Q b f + 0(a a ), (8) 

where N c — 3 is the number of colors. The experimental measurement, 8 

6R had (2m b ) = .36 ± .09 ± .03 (9) 

is in good agreement with the theoretical prediction from Q b em = — -|, veri- 
fying the SM electric charge assignment of the b quark. 

Similarly, the SU{2) L quantum numbers of the left- and right- handed 
b quarks are probed by the decay rate for Z to bb quark pairs. If the b did 
not have a top quark partner, it would be an isospin particle (T| = 0) 
and the decay width would be dramatically different from that of the SM. 
The decay width is given in terms of the left and right- handed couplings 
of the Z to the bb pair, 9 

r(Z^bb) = ^l(Ll + Rl). (10) 

If the b quark were an isospin singlet, then the decay width would be 
changed, 

T{Z ->bb) T *=- 1 _ 1 + 4Q fc sin 2 W + 8Q 2 sin 4 9 W 
T{Z^bb) T z=° ~ 8Q 2 b sm 4 9 w 
- 13 . 

The measurement of the Z — > bb decay width excludes the T| = hypoth- 
esis for the b quark. 3 

• The measured b couplings, combined with anomaly cancelation, 
require that the b quark have a T| = |, color triplet, fermion 
partner: this is the top quark. 
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Figurc2. Fcynman diagrams for gauge boson self-energies which give contributions pro- 
portional to M t 2 . 



2.3. Precision Measurements 

Before the top quark was discovered, an approximate value for its mass was 
known from precision measurements, which depend sensitively on the top 
quark mass. At tree level, all electroweak measurements depend on just 
three parameters: the SU{2)l x U(1)y gauge coupling constants and the 
Higgs vacuum expectation value, v — 246 GeV. These are typically traded 
for the precisely measured quantities, a, Gf, and Mz- All electroweak 
measurements at lowest order can be expressed in terms of these three 
parameters. Beyond the lowest order, electroweak quantities depend on 
the masses of the top quark and the Higgs boson. 

A typical example of the role of the top quark mass in precision mea- 
surements is the calculation of the p parameter, 



Ml cos 2 9 W 
Aww(O) A zz (0) 
Ml Ml ■ 



(11) 



where A vv is defined by the gauge boson 2- point functions, 

iU v v v = A vv (p 2 )g^ + B vv p»p» . (12) 

At tree level, the p parameter in the Standard Model is exactly one, but 
at one loop it receives contributions from gauge boson, Higgs boson, and 
fcrmion loops. The largest corrections are those involving the top quark 
loop. For simplicity, we compute only the corrections proportional to M 2 , 
which arc found from the diagrams of Fig. 2. 

The 2-point function of the Z boson is given by, 



d n k 

K 4c w ' vv J (27r)™ [den] 



m^ z = -N^fiif / —— (is) 
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Figure3. Fits to electroweak data. The bars show the pull from the Standard Model 
global fit. 3 



where cw = cos 9w and , 

T» v = Tr[(k + M t )-f(RtP+ + L t P_)(k -p+ M t )Y(RtP+ + L*P-)] 
den = [k 2 - M 2 ] [(it + pf - M 2 ] , (14) 

and P± — |(1 ± 75). The coupling of the Z to the top quark is given in 
Eq. 5 and p is the external gauge boson momentum. Shifting momentum 
in the integral, k' — > k — px, and using Feynman parameters, the integral 
in n = 4 — 2e dimensions becomes, 

iltt" = -N c (^-) 2 (i) 2 [ ^X- t dx 
zz ^c w ' v ; J (2tt)« J 

e 
2 



[4k' 2 (-l + § )(R 2 + L 2 ) + 16L t R t M 2 ]g^ + 
{k' 2 +p 2 x(l-x)-M 2 ) 2 



(15) 
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where we retain only those terms contributing to A zz (0)- The result pro- 
portional to M t 2 is, 



The analogous result for the W two-point function is, 



g 2 N c f 4tt , 1 



^-°>-s? U M 'U + 5i' (17) 



Combining Eqs. 16 and 17, the contribution to the p parameter which is 
proportional to M t 2 is found from Eq. II, 10 

g 2 N c M? G F N C M? 

^-64^M^-^7^' (W) 

Experimentally 4 p — 1.00126t;ooi4 4 an d so an upper limit on the top quark 
mass can be obtained from Eq. 18. 

Many of the precision measurements shown in Fig. 3 are sensitive to M t 2 
and so a prediction for the top quark mass can be extracted quite precisely 
by combining many measurements. In fact, precision measurements were 
sensitive to the top quark mass before top was discovered at Fermilab! The 
agreement between the direct measurement of the top quark mass in the 
Fermilab collider experiments and the indirect prediction from the precision 
measurements (as shown in Fig. 4) is one of the triumphs of the Standard 
Model. It is interesting to note that a similar limit on the Higgs boson mass 
from precision measurements gives Mh < 193 GeV at the 95% confidence 
level. 3 Precision measurements depend logarithmically on the Higgs mass, 
and so it is much more difficult to bound the Higgs mass in this manner than 
it is to restrict the top quark mass. Increases in the experimental precisions 
of the top quark mass and the W boson mass in Run II measurements at 
the Tevatron will provide an improved bound on the Higgs boson mass, 
SM h /M h ~ 40%. 11 

In the next sections, we discuss the discovery of the top quark at Fer- 
milab and the experimental exploration of the top quark properties, both 
at Fermilab and the LHC and at a future high energy e + e~ collider. 



3. Top Quark Properties 
3.1. Hadronic Production 

The top quark was discovered in 1995 at Fermilab in pp collisions at >/S = 
1.8 TeV. 1 This data set (called Run I) consists of an integrated luminosity of 
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Figure4. Measurements of the top quark mass at Fermilab (CDF and DO) and indirect 
predictions from precision measurements (LEP1, SLD and Myv)? 



L ~ 125 pb^ 1 . Both CDF and DO are currently rediscovering the top quark 
in the Run II data set. Run II will produce roughly 500 clean top quark 
events for each inverse femtobarn of data and so precision measurements of 
many top quark properties will be possible. 

In hadronic interactions, the top quark is produced by gluon fusion and 
by qq annihilation as shown in Fig. 5, 



99 

qq 



tt 
it. 



(19) 



The hadronic top quark production cross section, cry, at the Tevatron, 
pp — > it, (or pp — > it at the LHC) is found by convoluting the parton level 
cross section with the parton distribution functions (PDFs), 



<th(S) 



f l {xi,p)fj{x 2 ,n)a lj {x 1 x 2 S, fi), 



(20) 



where the hadronic center of mass energy is S, the partonic center of mass 
energy is s = XiX 2 S, the parton distribution functions are fi(x,p), and the 
parton level cross section is aij(s,fi). The parameter [i is an unphysical 
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Figure5. Fcynman diagrams contributing to top quark pair production in hadron col- 
liders. 

renormalization/factorization scale. To 0{a 2 8 ), the sum in Eq. 20 is over 
the qq and gg initial states shown in Fig. 5. 

It is useful to parameterize the parton level cross section (to next-to- 
leading order in a s ) as, 



Mi 



h%{p) + hlja a (p) + h%a B (ji) logf -^L 



Mi 



At lowest order, 0(a%), the parton level cross sections are 
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h°aa = , 



(2 + p) 



(3 



where 



(p 2 + 16/3 + 16) log 



AM} 



1 + /3 
1-/3 



-28-31p 



(21) 



(22) 



(23) 



The threshold condition is p = 1 and so &ij vanishes at threshold. 

To any given order in a s (/j,), the hadronic cross section gh must be 
independent of p, 

da H 



d log ps 



= 0. 



(24) 
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Applying this restriction to Eq. 21 yields, 







(df(xi,n) 

1 ~al 2"/(^2,m) + /(zi,m)- o, 2 

+ [ f(x 1 ,fj,)f(x 2 ,fJb) da ^,\ X2S o ^ dxxdx 2 . 



<t(xiX25, ii)dx\dx2 



aio gM 2 (25) 

The scale dependence of the PDFs is governed by the Altarelli-Parisi evo- 
lution functions, Pij{z): 



dfi(x,fJ.) = ^ 
9 log p 2 27r 15 



1 dy 



Pi. 



(26) 



Similarly the scale dependence of the strong coupling constant, a s (p), is 
governed by the evolution of the QCD (3 function. 

da s (p) 



d log p? 



-boat 



(27) 



where 6 = ^ 33 ^ 2 2 "' / ^ , and n// is the number of light flavors. At next- 
to-leading order the scale dependence of the various terms must cancel, 
yielding a prediction for ft?-, 
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(28) 



Using the explicit forms of the Altarelli-Parisi evolution functions, combined 
with Eq. 22, we find 12 
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(29) 
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where we have defined, 

M=Lh (^-L h (^. (30) 

The quantities hjj can only be obtained by performing a complete next- 
to-leading order calculation 12 ' 13 and analytic results are not available, al- 
though a numerical parameterization is quite accurate. The strong cou- 
pling evaluated at the top quark mass is small, a s (M t ) ~ O.f , and so a 
perturbative expansion converges rapidly. At energy scales significantly dif- 
ferent from the top quark mass, there are also large logarithms of the form 
(a s log(^i-))", which can be summed to all orders to obtain an improved 
prediction for the cross section. 1 The inclusion of next-to-leading- logarithm 
effects reduces the scale dependence of the cross section to roughly ±5%. 15 

Fig. 6 shows the Run I CDF and DO measured top production cross 
sections and top quark mass compared with two theoretical predictions. 
The Run I tt cross sections are, 16 

a tJ = 6.5±i;I pb CDF 

a a = 5.9 ±1.7 pb DO. (31) 

The curves labelled NLO+NLL include some of the logarithms of the form 
(a s log(^j-))" and the agreement with the experimental results is clearly 
improved from the NLO prediction alone. An updated theoretical study 
including next-to-leading-logarithm resummation gives the prediction 

^(VS = 1.8 TeV) = (4.81 - 5.29) pb for M t = 175 GeV, (32) 

with the range in the prediction corresponding to M t /2 < /j, < 2M t . 15 
The agreement between the predicted and the experimental production 
rates implies that the top quark is a color triplet, since that rate would be 
significantly different for a different color representation. 

At the LHC, the subprocess gg — ► tt dominates the production rate with 
90% of the total rate and the cross section is large, <jlhc ~ 800 pb, yielding 
about 10 s tt pairs/year. Because of the large rate, the LHC can search for 
new physics in top interactions. The large rate also has the implication that 
top production is often largest background to other new physics signals. 
A detailed understanding of the top quark signal is therefore a necessary 
ingredient of new physics searches. 
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Figure6. Fcrmilab Tcvatron Run I cross sections for pp — ► it. The dashed line is the 
NLO prediction, while the solid lines incorporate the next-to-leading logarithms. 14 



3.2. Weak Interactions of Top 

The weak interaction eigenstates are not mass eigenstates, 



L = - l5)V tq qW+ + h.c. (33) 

and so the interaction of the top quark with a W boson and a light quark is 
proportional to the CKM mixing element, V tq . Measurements of the decay 
rates of the top quark into the lighter quarks therefore translate directly 
into measurements of the CKM mixing angles. 

The unitarity of the CKM matrix and the restriction to 3 generations 
of fermions gives a limit on Vtb from the measured values of | V ub | and 
I V cb |: 4 

1 =| V ub | 2 + | V cb | 2 + | Va I 2 (34) 

giving, 

0.9991 <| V tb |< 0.9994 . (35) 

If there are more than 3 generations, however, Eq. 34 is no longer valid 
and there is almost no limit from unitarity on Vtb- 

Unitarity in the three generation model can be tested by measuring the 
ratio of the rate for top decays to the b quark to the rate for top decays to 
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Figurc7. Fcynman diagram for top decay to Wb. 



lighter quarks, 

_ T(t -^Wb) _ | V tb | 2 

tb ~ r(t - w q ) | v td | 2 + 1 y ts | 2 + | Vt 6 | 2 

= Mtf A (CDF) (36) 

This quantity can be measured at Fermilab by counting the number of 
tagged b's in a top quark event. Since the b quark lives approximately 
1.5 ps, the b quark will travel 450 am before decaying. The secondary 
b vertex can then be measured with a silicon vertex detector. Assuming 
unitarity of the CKM matrix, the denominator of Eq. 36 is 1, and the 
measurement can be interpreted as a measurement of Vtb, 

I |= .97t : & (37) 

consistent with Eq. 35. 

If there is a fourth generation of quarks, (t',b'), the charge — | b' is 
experimentally restricted to be heavier than m t — Mw, 4 and so the top 
quark cannot decay into the 6'. Then the denominator of R t b need not 
equal 1, and the measurement implies only | Vtb \»\ Vts |, I Vtd \- 

The direct measurement of Vtd m single top production will be discussed 
in Section 3.6. 



3.3. Top Quark Decay 

The top quark is the only quark which decays before it can form a bound 
state. This is due to the very short lifetime of the top quark, 

n~5x 1(T 25 sec , (38) 
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compared to the QCD time scale: 

t qcd ~ 3 x 1(T 24 sec . (39) 

This is very different from the b quark system, where the b quark combines 
with the lighter quarks to form B mesons, which then decay. 
Since | Vtb |~ 1, the dominant decay of the top quark is 

*(p)-> W»( Pw ) + b(p'), (40) 

The lowest order amplitude is shown in Fig. 7 and is, 

A (t - W"b) = --JL_u( p ')^(l - l5 ) u (p) . (41) 

Squaring the amplitude, and summing over the W polarization vectors, 

Y,e»{p w y»{p w ) = -<T + (42) 

gives the amplitude-squared, 

M 2 - 



2M 2 , 



(43) 



\A (t^Wb) \ 2 =g\M?-M^) 1 + 
The total decay width is 

r (t - wb) = Q)-^/ I Mt - wb) | 2 (dps 2 ) . (44) 

The factor of \ is the average over the initial top quark spin, and the phase 
space factor is, 

M 2 - M 2 

<™ = \& ■ (45 > 

Including the CKM mixing, the final result is then, 

ro(t ^ ) . |V .|.^( 1 -^)'(. + ^). <«, 

Higher order QCD corrections can be calculated in a straightforward man- 
ner and yield a precise prediction for the decay width, 9,18 

T(t -» bW+) =| V tb | 2 1-42 GeV. 
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Figurc8. Spin vectors for longitudinally polarized W's produced in the decay t — > Wb. 

3.4. W Helicity in Top Quark Decay 

The helicity structure of the top quark decays is interesting because it yields 
information on the V — A nature of the tbW vertex as is clear from Eq. 41. 
The top quark is produced through the weak interactions as a left-handed 
fermion (neglecting the b quark mass), which has spin opposite from the 
direction of the top quark motion. 

If the W is produced as a longitudinal W (helicity 0), then its momen- 
tum and polarization vectors are (see Fig. 8): 

Pw = (£V,0,0, | Pw I) 

ef = j^(\pw\,0,0,E w ) 



\Jk1 + (^-). (47) 
M W \ M t J 1 ' 



The amplitude for a top decaying into a longitudinal W is then, 
A(t _ bW L ) = --^=5^(1 - 75)t + (Mw 



2^2 M w y '°' V M t 

This gives the result that the decay of the top into longitudinal W's is, 

T(t -> bW L ) ~ ^ . (49) 

The decay of the top quark into a positive helicity W, h w = +, is 
forbidden by angular momentum conservation since a massless b quark is 
always left-handed, while the heavy top can be either left- or right-handed. 
This is illustrated in Fig. 9. CDF has measured the decay of the top quark 
to a right handed W and finds a result consistent with 0, 19 

BR(t^> bW R ) = .11 ±0.15 . (50) 

For a transversely polarized W, the polarization vector is given by, 

e^ = i=(0,l,±i,0). (51) 
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Figurc9. Spin vectors for positive helicity (right handed) Ws produced in the decay 
t — > Wb. Conservation of angular momentum forbids this decay. 
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FigurclO. Spin vectors for transverse helicity W's produced in the decay t — » Wb. 



This gives a decay rate, 

T(t -» bW T ) ~ g 2 m t . (52) 

Comparing Eqs. 49 and 52, it is clear that the rate for top decay to trans- 
verse W's is suppressed relative to that for decay to longitudinal W's, 

r(t - 6w L ) 

-f'o = 



1 



2M^ 

= 70.1% . (53) 

Can we measure this? The W helicity is correlated with the momentum 
of the decay leptons: h w = + gives harder charged leptons than h w = — . 
A preliminary measurement of the polarization of the W's in top decay 
yields, 19 

F = 91 ± 37 ± 13% CDF, (54) 

consistent with Eq. 53. With 2 fb^ 1 of data from Run II at the Tevatron, 
a more definitive (~ 5%) measurement of Fq will be possible. 20 
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3.5. Top Production and Decay 

The Tevatron and the LHC produce predominantly it pairs, which decay 
almost entirely by t — > Wb (since | V t b |~ l), a 

pp^tt^W + W-bb. (55) 

It is useful to classify top quark events by the decay pattern of the W + W~ 
pair: 

• Di-lepton: Roughly 5% of the it events fall into this category, 

W + ->l + v, l = e,u 

W~ -» l~v . (56) 

Since Br(W -» ev) + BR(W -» ^) ~ 22%, the 125 p^ 1 at 
Run I of the Tevatron produced <~ 625 it events which yielded 
<~ .22 x .22 x 625 <~ 30 di-lcpton events. These events are clean, 
but suffer from the small statistics. 

• lepton + jets: This class of events comprises roughly 30% of the 
total, 

W ± -» l±v, W T -» jets . (57) 

Lepton plus jet events are fully reconstructable and have small 
backgrounds. Assuming at least 1 b quark is tagged, CDF ob- 
served 34 such events in Run I, with an expected background of 
8. 22 In Run II, with 2 fb' 1 , there should be 1000 lepton+jets 
events/experiment. This increase from Run I results from higher 
luminosity, the higher cross section at VS = 2 TeV compared to 
=1.8 TeV, and upgraded detectors with improved capabilities 
for ^-tagging. 

This channel produced the most precise measurement of the 
cross section for top pair production. The combined average from 
CDF and DO at Run I is, 

a(pp^ti) \^s =1 . STeV = 5.9 ± 1.7 pb (58) 

• All jets: 44% of the events fall into this category. This class of 
events was observed by both CDF/DO at Run I, with large back- 
grounds. 



a Rcccnt reviews of the experimental issues involved in top quark physics can be found 
in Ref. 21. 
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Table2. Contributions to single top 
production at the Tevatron and the 
LHC. 





Tevatron, Runll 


LHC 




2.1 pb 


238 pb 


O-qq 


.9 pb 


10 pb 


O-gb 


.1 pb 


56 pb 




Figurcll. Feynman diagrams contributing to single top production at the Tevatron and 
the LHC. 



3.6. Single Top Production 

Single top production can provide a precise measurement of Vtb- The dom- 
inant production processes at the Tevatron are the t channel W exchange 
and the s channel qq annihilation shown in Fig. 11. At the LHC, the sub- 
process gb — > tW is also important. The contributions (to NLO QCD) of 
the various subprocesses at the Tevatron and the LHC are given in Table 
2 23,24 f these processes are proportional to | Vtb \ 2 ■ 

Single top production was not observed in Run I, and both experiments 
set limits on the production rate. At Run II of the Tevatron, and at the 
LHC, single top production can be observed with a rate roughly half that 
of the it rate and will lead to a measurement at the Tevatron of, 

I SVtb I W 
I — |~5% 

with 30 fb^ 1 of data. With 2 the Tevatron Run II should achieve an 

accuracy of | |~ 13%. 20 

3.7. Measurements of Mt 

M t is a fundamental parameter of the Standard Model, but the value is 
not predicted (except indirectly from precision measurements such as the 
p parameter discussed earlier). In the Standard Model, a precise value 
of M t is important primarily for its role in precision electroweak observ- 
ables and for the prediction for the Higgs boson mass extracted from these 
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Figurel2. Implications of precision measurements of Mt and M\y on the indirect ex- 
traction of the Higgs boson mass. The Standard Model allows only the shaded region. 3 



measurements. In extentions of the Standard Model, such as those dis- 
cussed in Section 5, a precise knowledge of M t plays a critical role in defining 
the parameters of the theory. 

The top quark mass was measured in Run I at the Tevatron and the 
most accurate value of M t comes from the single lepton+jets channel, 2 

Mt = 176.1 ± 4.8 (stat) ± 5.3 (syst) CDF 

= 173.3 ± 5.6 (stat) ± 5.5 {syst) DO . (59) 

In the lepton plus jets channel there is one unknown parameter (the longitu- 
dinal momentum of the v coming from the W decay) and three constraints 
from the reconstruction of the W masses and the requirement that the re- 
constructed masses of the top and anti-top be identical, M t = Mj. The 
best value of the top quark mass found from combining all channels at the 



21 



Tevatron is, 

M t = 174.3 ±5.1 GeV . (60) 

The systematic error is dominated by the uncertainty on the jet energy 
scale, 

±4.4 GeV CDF 
±4 GeV DO . (61) 

With increased data from Run II, the mass measurement at the Tevatron 
will be improved to SM t ~ 3 GeV with an integrated luminosity of 2 

Fig. 12 shows the impact of a precise measurement of the top quark 
mass combined with a W mass measurement for the prediction of the Higgs 
boson mass. Clearly, for the precision measurements of the W and top 
masses to be useful, the experimental error on both measurements needs to 
be comparable. The Higgs mass prediction extracted from precision mea- 
surements can be compared with the direct measurement of the Higgs mass 
at the LHC and the consistency of the Standard Model tested. Tevatron 
Run II with 2/6 -1 should achieve accuracies: 

SM W ~ 27 MeV 
SM t ~ 3 GeV , (62) 

yielding a prediction for the Higgs mass with an uncertainty of 25 

^ ~ 40% . (63) 

The LHC will have large statistics on top quark production due to the large 
rate and can use the lepton+jets channel with 10 fb^ 1 of data to obtain a 
measurement 26 

SM t ~ .07 GeV . (64) 



3.8. Top Spin Correlations 

When pair produced in either e + e~ or qq collisions, the spins of the ti 
pair are almost 100% correlated. 27 The top decays before the spin can flip, 
and so spin correlations between the t and i result in angular correlations 
among the decay products. At threshold, the ti cross section is s-wave, and 
the e+e~ or qq spins are translated directly to the ti spins. The incoming 
q,q have opposite helicities, since they are interacting through a gauge 
interaction and the ti spins aligned along the beamline. At high energy, 




y/s >> m t , the top mass is irrelevant and the t,i must have opposite 
helicities. Finding the correct basis for intermediate energies is subtle, 
however. 



4. The Top Quark as a QCD Laboratory 

The production and decay mechanisms for the top quark in both e+e~ and 
hadronic collisions are well understood and so the top quark system pro- 
vides an ideal laboratory for studying pcrturbative QCD. This section con- 
tains a pedagogical treatment of the 0(a s ) next- to- leading order (NLO) 
QCD corrections to the process e+e~ — > it. 28 This process is particularly 
simple since the QCD corrections affect only the final state. The calcula- 
tion is performed using the two cut-off phase space slicing (PSS) algorithm 
for evaluating the real gluon emission diagrams. 29 ' 30 The total rate for 
e + e~ — > it tests our understanding of pcrturbative QCD and additionally 
the threshold behavior of the total cross section can be used to obtain a 
precise measurement of M t , along with information on a s (M t ) and the tih 
Yukawa coupling. An understanding of the QCD corrections are crucial for 
these interpretations. 

As a second example of the role of QCD in top quark physics, we con- 
sider the associated production of a it pair along with a Higgs boson, h. 
The tih Yukawa coupling is a fundamental parameter of the theory and 
probes the mechanism of mass generation, since in the Standard Model, 
gtth — -^7- The processes e + e~ — ► tih and pp(pp) — ► tih provide sensitive 
measurements of the top quark- Higgs coupling and can be used to verify 
the correctness of the Standard Model and to search for new physics in the 
top quark and Higgs sectors. The importance of the QCD NLO corrections 
to tih production is discussed at the end of this section. 
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4.1. NLO QCD corrections to e+e — > tt 

The Born amplitude for e + e~ — > proceeds via s— channel 7 and Z ex- 
change, as shown in Fig. 13. The dominant contribution at y/s = 500 GeV 
is from 7 exchange and so for simplicity we will neglect the Z contribution 
here. (The complete 0(a s ) result, including the contribution from Z— ex- 
change can be found in Refs. 28 and 30.) The electron mass can also be 
neglected. The lowest order amplitude is: 

A LO (e + e- tt) = e ^^lv{ q )Yu(p)u(p')^v( q ')8 l3 , 

where i,j are the color indices of the outgoing top quarks. It is simplest to 
work in the center-of-mass frame of the incoming e + e~ pair, where the 4— 
momenta of the particles are given by: 

p= ^(1,0,0,1) 
q= ^(1,0,0,-1) 



p' = ^(l,O,/3sin0,/3cos0) 
q' = ^(1, 0, -/? sin 6,-13 cos 6) 



with 



S = (p + q) 



2 



, 3 = V 1-M (65) 

and 9 is the angle of the top quark emission with respect to the electron 
beam direction. The threshold condition is defined by j3 — 0. 

It is straightforward to find the tree level amplitude-squared in n = 4— 2e 
dimensions, (summed over colors), 



A LO | 2 = 64ir 2 a 2 Q 2 t N c 



2(l-e)+/3 2 (cos 2 0- 1) 



(66) 



The color factor is N c = 3. The corresponding total cross section is: 

aLO = h J 1 Alo |2 (dP5)2 - (67) 

Integrating over the two body phase space, 

o2e /a \e 1 fir 

^ = wA-)^w^)L sin1 "^' (68) 
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*(«') 



Figurel4. Vertex correction to e+e — > it. The 7*" denotes an off-shell photon which 
is connected to the initial e + e~ state. 

and including a factor of \ for the average over the initial electron spins, 
gives a total rate, 



The 0(e) terms which we have neglected will be important in computing 
the 0(a s ) QCD corrections. The cross section vanishes at threshold, (3 — > 
0, and decreases with increasing energy, a <~ -. This has the obvious 
implication that the higher energy an e + e~ collider, the larger the required 
luminosity. 

4.2. 0(a s ) corrections to e+e - — > it 

The 0(a s ) corrections to e+e~ — > ti contain both real and virtual contri- 



The one-loop virtual contribution consists of a vertex correction and a wave- 
function renormalization of the top quark leg. The vertex correction is 
shown in Fig. 14 and can be parameterized as: 




2irN c a 2 Q 2 t 
3s 



/3(3-/3 2 ) + 0(e). 



butions: 



&NLO = &LO + O virtual + CreaZ- 



(69) 



A 



vertex 



= eQ e v{q)lvu(p)-T»(p',q') 



(70) 



with. 



r„(pW) = (-eQ t )6, 



l3 glC F [N\u(p') A 




+ B lv v(q') 



(71) 
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where the normalization is, 

[N] = 

and 



16tt 2 



47T/X 2 



C F = 



2N r 



r(l + e) 



(72) 



(73) 



is the color factor for the vertex diagram. The parameter fj, is an arbitrary 
renormalization scale. Explicit calculation gives the results, 28 ' 30 

1 



A = {fi- 
A/3 = log 



1-/3 
1 + /3 



5=1 
e 



313 A/3 + 



4 A ^ + 2A /3 l0 §(j 



2/? 



+ /3 



+2Li 2 ( 



1-/3, 2tt 2 
-) H 

i + 3 



= — + -D/m • 

The di-logarithm is defined by, 
Li 2 (x) 



dz- 



1 — xz 



(74) 



The 0(a s ) contribution to the virtual cross section is the interference 
of the vertex correction with the lowest order result. (The A term of Eq. 71 
does not contribute to the interference with the lowest order result when 
the electron mass is neglected). 



^vertex 



«s / 47T^ 

27 V~M? 



(dPS) 2 A LO A* vertex 
C F T(l + e)a LO 



B * + R 



(75) 



Note that the (e) terms from the lowest order amplitude and from the 2- 
particle phase space combine with the term to give a finite contribution 

to Cyertex • 

The wavefunction renormalization for the top 
quark must also be included in the virtual correction, as shown in Fig. 
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tip') 




tiq') 

Figurcl5. Fcynman diagram for top quark wavcfunction rcnormalization. 



15. This diagram gives a multiplicative correction to the lowest order rate. 
Using the on-shcll rcnormalization scheme, 



az = 25Z tt <JLO, 



(76) 



where 



dZu = Q __ 9 | P 



dp 2 



4ttV M 2 



Mr 



C F r(l + e) 



The total virtual contribution is thus: 28 ' 30 



(77) 



&VIRT = O vertex + (T Z 

= a L0 ^-T{l + e)C F 

ZlT 



\M 2 



{-3 + B e )~ -4 + B/i 
e 



(78) 



4.3. Real contributions 

The virtual singularities of Eq. 78 are cancelled by those arising from real 
gluon emission, e+e~ — ► tig. A sample diagram of this type is shown in 
Fig. 16. The real gluon emission diagram of Fig. 16 contains a fermion 
propagator on the outgoing top quark leg of the form: 
1 1 1 

[p 1 +P g ) 2 -M 2 ~ 2p> ■ p g ~ 2E'E g (l-0cos§y 

where E 1 and E g are the energies of the outgoing top quark and gluon, 9 is 
the angle between them, and p = y 1 — -^rj . Because of the massive top 
quark, (3^1, and this diagram has no collincar singularity as 9 — > 0. There 
is, however, a soft singularity as the energy of the gluon vanishes, E g — ► 0. 
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The soft singularity is regulated using a technique called phase space 
slicing (PSS). 30 This method divides the gluon phase space into a soft gluon 
plus a hard gluon region: 

&real — &soft ®hard • 

(79) 

The soft region is defined by the energy of the emitted gluon, 

E g <S s ^-, (80) 
while the hard region is the remainder of the gluon phase space. 

E 9 > • (81) 

The separation into hard and soft regions depends on the arbitrary param- 
eter S s . Individually, a so f t and <7hard depend on the cutoff, 5 Sl and the 
independence of the sum is a check on the calculation. 5 S must be small 
enough that terms of 0(S S ) can be neglected, while large enough to pre- 
vent numerical instabilities. Ohard is finite and can be evaluated numerically 
using standard Monte Carlo techniques. 

In the soft gluon regime, the cross section can be evaluated analytically. 
The three body phase space is evaluated in the soft limit in the rest frame 
of the incoming e + e~: 

p + q = p' + q' + Pg , (82) 

giving a gluon energy of 

Spi q i 



2^7s 
{p' + q'f ■ 



(83) 
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The 3-body phase space in the soft limit can be separated into a 2 body 
phase space factor times a soft phase space factor: 31 

f d n ~ x r>' I" d"~ 1 o' 

{dFSh = J / ■ l2 ' ) " 6 " (p + - 

(84) 



/ 



2E g (2n)"-i ■ 



For the soft limit, we set the gluon momenta p g = in the delta function 
and the phase space then factorizes: 

(dPSHsoft) = (dPS) 2 J ^r^zi • (85) 

The phase space is most easily evaluated by choosing an explicit repre- 
sentation for the gluon momentum, 

p g = E g (l, ...,sin#i sin 2 , sin#i cos#2, cos#i). (86) 

yielding, 32 

d n - 1 Pg =d\p g \\Pg I"" 2 dn n _ 2 

= dE g E n g - 2 sin™" 3 BydQx sin™- 4 9 2 d6 2 dn n - 4 (87) 

where, 

27r (n-3)/2 

^n-A = • ( 88 ) 

The integrals can be performed analytically to find the three body phase 
space in soft limit: 32 

(dPSHsoft) = (dPS) a (±E )\^)§^dPS soft 



1 Is 



e rS s 



^2 



dPS *oft = ^) J o dE g E]^ 

■ sin 1_2e 9-idB-i. sin _2e 9 2 d8 2 . (89) 

Note that the soft phase space depends explicitely on the cut-off S s through 
the limit on the E g integral. 
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The soft phase space must be combined with the amplitude for e+e 
tig in p g — > limit: 



isoft 



e 2 Q e Q 



-9 s T A v{q)-i^u{p) 
s 



g s T J 



P' 'Pg l' ■ Pg 



Alo 



(90) 



where we have set p g = in the numerator. The approximation of retaining 
only the most singular contributions is known as the eikonal ( or double 
pole) approximation. 33 The soft limit of the amplitude always factorizes in 
this manner into a factor multiplying the lowest order amplitude. 
The soft amplitude-squared is 



A soft \ 2 = g 2 Tr(T T ) \ A LO 



M 2 



M 2 | s - 2M 2 



(p' -PgY (<?' -Pg) 2 ' P' -PgQ' 'Pg 

giving the cross section in the soft limit, b 

Usoft = h{\) J 1 Asoft |2 ( dPS ^ s °w ■ 



(91) 



(92) 



The problem reduces to evaluating the soft integrals. A typical integral is 

1 



jsoft 



(P' ■ Pg) 



;dPS 3 ( S oft) . 



(93) 



Since the phase space and the amplitude squared are Lorentz invariant, we 
are free to chose any convenient frame. We use 



p g = E g (l, ... sin 0i sin # 2 , sin#i cos# 2 , cos#i) 



0,(3) 



(94) 



b As always, the factor of 1/4 is the spin average. 
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and the soft integral is, 



7T V 4 



dEgE 1 - 2 '- sin 1 " 26 6 l d6 l sin" 2 * d0 2 



\sE*J (1 -f3cos6 1 ) 2 



2e 



(95) 



The complete set of soft integrals relevant for this process is given in the 
appendix of Rcf. 34. 

The soft cross section is 

'4VY r (!- e ) 



a soft = olo-z-Cf 
Ztt 



T(l - 2e) 



O e 

r W* 

e 



(96) 



where B e — —C e so the sum of virtual plus soft contributions is finite. The 
finite contribution, C/j n , depends on S s and is given by 30 



Cfi 



-4log(S s )--A -2[-+P 



Li 2 



2/3 
1 + /? 



A 2 



/3 

■ it is then, 



0(5.) 



(97) 



The complete 0(a s ) result for e+e 

ffJVLO = ^LO + ^virtual + soft + &hard ■ (98) 

The combination <J S oft+&hard is independent of S s , while the | poles cancel 
between cr so /t + ^virtual- 

This example, e + e~ — » tt, is simple enough that it could have been 
done analytically without the approximation of phase space slicing 28 . Phase 
space slicing, however, is useful for more complicated examples and also to 
construct Monte Carlo distributions beyond the leading order. 



4.4. Threshold Scan in e+e — > tt 

The energy dependence of the cross section for e+e~ — > it at threshold 
depends sensitively on the parameters of the top sector. The location of 
the peak depends on M t , while the height of the peak depends on T t . The 
problem is that the location of the peak as a function of center-of-mass 
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Figurel7. Cross section for e + e — — > tt as a function of the center-of-mass energy in 
Ge V using the pole mass definition for the top quark mass at lowest order ( dot-dashed) , 
NLO (dashed), and NNLO (solid). The vertical axis shows <r(e+e — — » tt)/a p t, where 
Opt = Ana 2 /3s. The three sets of curves at each order represent the variation with the 
renormalization scale fi, (fi = 15, 30, and 60 GeV). 37 



energy shifts at higher orders and there is a large renormalization scale 
dependence to the prediction as seen in Fig. 17. 35 

At threshold, j3 — ► 0, the theory possesses two small parameters: a s 
and (3 and the cross section diverges at 0(a s ), 9 

a(e + e- ^tt)\ ^ o =a LO (l + ?^y (99) 

This divergence as f3 — > is known as a Coulomb singularity and non- 
relativistic quantum mechanics can be used to sum the leading contribu- 
tions in powers of ^. 36 

Beyond lowest order in a s , the predictions are quite sensitive to the pre- 
cise definition of the top quark mass. The most straightforward definition 
is the pole mass appearing in the top quark propagator, 

D & = ITT ?n ■ ( 10 °) 

P - M t [pole 

The pole mass, M t | po /e, is measured by reconstructing the four-momenta 
of top quark decay products. This definition of the top quark mass is, 
however, uncertain by QCD hadronization effects, O(Kqqd)- These effects 
connect the top quark with its decay products. Even more disturbing is 
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Figurel8. Cross section for e + e — — > tt as a function of the center-of-mass energy in 
GeV using the IS mass definition for the top quark mass at lowest order ( dot-dashed), 
NLO (dashed), and NNLO (solid). The vertical axis shows <r(e+e — — > tt)/a p t, where 
Opt = 47ra 2 /3s. The three sets of curves at each order represent the variation with the 
renormalization scale fi, (fi = 15, 30, and 60 GeV). 37 



the fact that the location and height of the peak in the e + e~ — > tt cross 
section change significantly between lowest order, NLO, and NNLO as seen 
in Fig. 17. A better definition of the top quark mass is the short distance 
IS mass: 37 ' 38 



mis = M t \pole (1 - 




(101) 



Using this mass definition, the cross section for e + e — > tt peaks at 2mis 
and the location of the peak does not shift at NLO or NNLO, as seen in 
Fig. 18. 

The NNLO result still has roughly a 20% scale uncertainty. It is esti- 
mated that by scanning the threshold dependence of the cross section for 
e + e~ — > tt a measurement with the accuracy 8M t ~ 200 MeV can be made 
with an integrated luminosity of L = 50 fb^ 1 . This is to be compared with 
the expected accuracy of the LHC measurement, SM t ^1—2 GeV. 26 A 
precise measurement of the top quark mass is of interest primarily because 
of its implications for electroweak precision measurements. 25 
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4.5. tth production in e+e and pp collisions 

In many models with physics beyond the Standard Model, the Higgs cou- 
pling to the top quark is significantly different from that in the Standard 
Model. In the Standard Model, the fermion-Higgs Yukawa couplings are, 

L=-g ffh W>h (102) 

where 

M f 

9ffh = -f ■ (103) 

Since the top quark is the heaviest quark, it has the largest Yukawa coupling 
to the Higgs boson. It is important to measure this coupling accurately in 
order to verify that the Higgs mechanism is the source of fermion masses. 

The Higgs boson couplings to the lighter quarks can be tested by mea- 
suring the Higgs decays to fermion pairs. This is not possible, however, for 
the top-fermion coupling. Only for very heavy Higgs bosons {Mh > 2M t ) 
is the decay h — > it kinematically accessible, and the branching ratio of the 
Higgs into it pairs is significantly smaller than that into gauge boson pairs. 
It does not appear feasible to measure g tt h through this channel. The most 
promising prospect for measuring the top quark Yukawa coupling is tth 
production in e+e~ and pp (or pp) collisions. 

4.6. tth at an e+e — collider 

At an e + e~ collider, tth associated production occurs through s-channcl 7 
and Z exchange, as shown in Fig. 19. The contribution from Z exchange 
is a few percent of the total rate and will be neglected here. With this as- 
sumption, the rate is directly proportional to g\ th - In order to interpret the 
cross section as a measurement of gtth, however, it is necessary to include 
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the QCD corrections which are potentially of the same order of magni- 
tude as new physics effects which may change gtth from the SM prediction. 
The NLO QCD result includes virtual and real gluon corrections on the 
outgoing t,t legs. The calculation follows that of the previous section for 
e + e~ — > tt, although in this case the virtual calculation also includes box 
diagrams. These corrections have been found in Refs. 39 and 40 and can 
be parameterized as 

K=^^ = l + a s {,i)F{s,MlMl) . (104) 

The parameter /i is an unphysical renormalization scale and, at this order, 
the only dependence on [i is in the running of a s ((i). For Mh = 120 GeV, 
and \x ranging from M t to y/s, 

K(^=500 GeV) ~ 1.4-1.5 

K(y^ = 1 TeV) - .8 - .9 . (105) 

The NLO results for tih production in e + e~ collisions are shown in Fig. 20. 
The rate is quite small and has a maximum value for y/s <~ 700 — 800 GeV 
for Mh in the 120 GeV region. 

In order to measure the top quark Yukawa coupling in e + e~ collisions, 
we consider the final state 

e+e- -» tih -» W + W-bbbb . (106) 

The events are then classified according to the W decays, in an analo- 
gous manner to that discussed in Section 3.5 for top quark decays. Since 
the event rate is quite small, it is advantageous to combine hadronic and 
semi-leptonic W decay channels. 41,42 At yfs = 800 GeV , an integrated 
luminosity of L = 1000 /6 _1 , can measure the coupling to an accuracy of 41 

^ ~ 5.5% . (107) 
gtth 

At a lower energy, yfs = 500 GeV, the rate is considerably smaller and 
even with an integrated luminosity of L = 1000 /6 _1 , the precision is quite 

42 

poor, 

^ - 20% , (108) 

gtth 

for M h = 120 GeV. 
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e + e -> tth 




M h =110 GeV 
M„=130 GeV 
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Figure20. NLO rate for e+e — + tth as a function of center of mass energy.' 



tth (or pp — ■+ tth) can be calculated in a 
tth. In the hadronic case, however, 



4.7. pp -> tth at NLO 

The QCD corrections to pp 
similar fashion to those for e + e" 
the calculation is complicated by the existence of strong interactions in 
the initial state. In this section, we discuss the treatment of initial state 
singularities in hadron collisions. The ingredients of an NLO calculation 
for pp or pp — > tth are: 34,43 

• ^virtual ■ This contains all of the one-loop diagrams contributing 
to the process. Finite diagrams, or diagrams which contain only 
ultraviolet singularities, can be evaluated numerically. 44 Diagrams 
which contain infrared singularities must be evaluated analytically. 



36 
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Figure21. NLO QCD corrections to pp -> tth at y/S = 2 TeV. This figure shows the 
independence of the total cross section on the soft cutoff, <5 S . 34 

• &soft'- This contains the soft singularities arising from both the 
initial and final state soft gluon radiation in the limit p g — > 0. The 
calculation is analogous to the discussion of Section 4.3. 

• Chard '■ This is the contribution from real gluon radiation when the 
gluon has an energy above the soft cut-off. In contrast to the e + e~ 
example of Section 4.3, this contribution has a singularity when 
the initial state quark or gluon is parallel to the radiated final 
state gluon. The initial state collinear singularities are absorbed 
into the parton distribution functions (PDFs). 

The combination a so f t + cr hard + a virt uai is finite. 

4.8. Collinear Singularities and Phase Space Slicing 

To illustrate the effects of collinear singularities, we consider the parton 
level subprocess, 

q(p) + q(q) ^ t(pt) + t{Pf) + h( Ph ) + g ( Pg ) . (109) 

As in the e + e~ example, we divide the gluon phase space into a soft and 
a hard region. We then further divide Ohard into a region with collinear 




divergences and a finite region, 

Ohard = & hard \coll +&hard \not-coll (HO) 

The separation depends on a new cut-off, 5 C , although a hard must be inde- 
pendent of 5 C . If the emitted gluon is from the initial state quark line, then 
there is a quark propagator, 

11 1 . . 

(Ill) 



(P-Pg) 2 ^P-Pg sfsE g (l - cos6) 
which becomes singular when p and p g are parallel (cos 9 — > 1). The 
hard/not-collinear region is defined by, 

m> s < 

In this region the quark propagator of Eq. 1 1 1 is finite and the cross section 
can be computed numerically. 
In the collinear region, 

2 P-Pg ^ x 2 1'Pg ^ x , 1lQ v 

E g V s E g V s 

and the matrix element factorizes into a parton i splitting into a hard parton 
i' plus a collinear gluon with kinematics 

Pi' = zpi 

Pg = (1 - z )Pt 

Sig = 2>Pi ' Pg 

These kinematics are illustrated in Fig. 22. In the current example, q — > 
and the matrix element factorizes, 34 

2P--,(z) 

| -» tt/i + g) \ 2 ^ co mnear 5V<£ I ^L («'j ~» «&) P " (H4) 

ZSia 
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where Pu> is the probability that parton i splits into parton i', with fraction 
z of the initial momentum. In this example, 

v2 



P iil {z) = P qq {z) = C F 



1 



1 



6(1-3) 



(115) 



The phase space also factorizes in the collinear limit: 31 



{PS)i{ij -> tt/l + g) collinear 

■ z dz ds ig [(l - z)s ig ] 



(PS) 3 (i'j -> ttfc) 



r(l - 2e) 16tt 2 



(1-z) 



Pi' • PjS c - Si 



(116) 



Combining Eqs. 115 and 116, the hard-collinear parton level cross section 
becomes: 45 



= a s r(l-e) 1 u 



j)j . 



(1 - z) 2 p 4 / -Pj 



zM 2 



Pu>(z)<JLo(i'j — > tt/i) 



(117) 



This result has a dependence on both <5 C and (5 S . The 1 — 5 S limit on the 
z integration excludes soft gluon region. The remaining singularities in 
Eq. 117 are absorbed into the mass factorization as described in the next 
section. 



4.9. Mass Factorization 

The hadronic cross section depends on the PDFs, fi(x): 



C7 H (S) = £ y - J f i (xi)fj(x 2 )cTij(xiX2S)dx 1 dx2 , 

where a is the parton level cross section. At lowest order, there is no 
dependence on the renormalization/ factorization scale, fi. Since the only 
physical quantity is the hadronic cross section, <jh(S), we can define scale 
dependent PDFs to absorb the hard/collinear singularity of Eq. 117, 30 ' 31 



fi{x,n) = fi(x) + 



a 



l\a s r(l-e) /4tt^ 
~e) 2^r(l-2e) \~t^~ 



1 ~ S " Ay T Q 1 

-PiAz)f v {-) + C F {2\og{5 s ) + -)\ . (118) 
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The upper limit on the z integration excludes the soft region. After con- 
voluting the parton cross section with the renormalized quark distribution 
function of Eq. (118), the infrared singular counterterm of Eq. (118) exactly 
cancels the remaining infrared poles of cr„j r t + o'soft + <?ha,rd/coll- 

The cancellation of the cutoff dependence at the level of the total NLO 
cross section is a very delicate issue, since it involves both analytical and 
numerical contributions. It is crucial to study the behavior of (Jnlo m a 
region where the cutoff(s) are small enough to justify the approximations 
used in the analytical calculation of the infra-red divergent soft and collinear 
limits, but not so small to give origin to numerical instabilities. This effect 
is illustrated in Fig. 21. 

The NLO result for pp — > tih at the Tevatron is shown in Fig. 23 
as a function of the unphysical factorization/ renormalization scale. The 
complete calculation includes not only the qq initial state described here, 
but also the gg and qg initial states. Note the reduced /i dependence of the 
NLO rate from the lowest order prediction. The cross section is very small, 
making this a challenging measurement. 

At the LHC, the process pp — > iih has a significant rate for Mh < 
130 GeV and we can look for h — > bb, giving the final state of W + W~bbbb. 
This is an important discovery channel for a Higgs with Mh < 130 GeV. 
The ATLAS experiment, with L = 100 fb' 1 , and assuming M h = 120 GeV 
can measure 5gtth <~ 16%. This is illustrated in Fig. 24. 26 The signal and 

9tth ° ° 

the background have similar shapes at the LHC, increasing the difficulty of 
this measurement. 



5. Heavy Top as Inspiration for Model Builders 

We turn now to a discussion of the top quark in models where there is new 
physics involving the top quark beyond that of the Standard Model. The 
fact that the top quark is much heavier than the other quarks gives it a 
special role in electroweak symmetry breaking. In this section, we examine 
the role of the top quark in the minimal supersymmetric model and in 
models with dynamical symmetry breaking. We pay particular attention 
to how data from the Tevatron and the LHC can help us to test these ideas. 

The large value of the top quark mass has important consequences in 
a supersymmetric (SUSY) model. 46 In the SM, the Higgs boson mass is a 
free parameter, while in a SUSY model there is an upper limit on the Higgs 
mass. At tree level in the minimal SUSY model, Mh < Mz, but the large 
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Figure23. Leading order and next-to-leading order rate for pp — * tth at the Teva- 
tron for M h = 120 GeV and \/S = 2 TeV as a function of the unphysical factoriza- 
tion/renormalization scale. 34 



top quark mass, and the subsequent large coupling to the Higgs boson, 

rru v2 

— 1, (119) 

v 

lead to large radiative corrections to the prediction for the Higgs boson 
mass, raising the lower limit on the Higgs mass significantly. In addition, 
corrections due to the large top quark mass typically imply that the scalar 
partner of the top, the stop, is the lightest squark. 

A further difference between the SM and a SUSY model which is due 
to the large top quark mass is the understanding of the mechanism of elec- 
troweak symmetry breaking in a SUSY model. If the supersymmetric model 
is embedded in a grand unified model, then the simplest version (mSUGRA) 
has the interesting property that the evolution of the parameters from the 
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Figurc24. Invariant mass for the bb pair from the Higgs decay in the process pp — > tth 
at the LHC with Mj, = 120 GeV (data points). The dashed line is the background. 
The error bars correspond to the statistical uncertainty with an integrated luminosity of 
L = 100 /fc" 1 . 26 

GUT scale to the weak scale drives the Higgs mass-squared negative as a 
consequence of the large value of M t , and therefore explains electroweak 
symmetry breaking. In fact, this mechanism requires that the top quark 
have a mass near its observed value. 

In models with dynamical symmetry breaking, the large value of the 
top quark mass can imply that the top quark is different from the lighter 
quarks, perhaps having different gauge interactions in an extended gauge 
sector. In Section 5.6, we discuss topcolor as an example of a model of dy- 
namical symmetry breaking where the top quark plays a special role. Many 
possibilities for dynamical symmetry breaking in electroweak interactions 
are detailed in the review by Hill and Simmons. 47 

5.1. EWSB and the top quark mass 

We begin by considering the importance of the top quark mass in the SM. 
In the SM, electroweak symmetry breaking occurs through the interactions 
of an SU{2)l doublet scalar field, with a scalar potential, 48 
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V = /j, 2 $& + A($$t) 2 . (120) 

If /Lt 2 < 0, the neutral component of the Higgs doublet gets a vacuum 
expectation value, 

(*0> - ^ (121) 

and the scalar potential has the shape shown in Fig. 25. Electroweak 
symmetry breaking (EWSB) occurs and masses are generated for the W 
and Z gauge bosons. After EWSB, there remains in the spectrum a physical 
particle, the Higgs boson, h, and at tree level the quartic self- coupling of 
the Higgs boson, A, is related to the Higgs mass by 

Ml = 2Xv 2 . (122) 

A measurement of this relationship would test the mechanism of EWSB. 

Eq. 120 is the scalar potential at the electroweak scale, A <~ v. All 
parameters of the electroweak theory change with the energy scale A in a 
precisely calculable way, and so will be different at higher energies. The 
most interesting parameters for us are the Higgs quartic self-coupling, A, 
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and the top quark Yukawa coupling, g tt h — 4^ • To one loop accuracy, 



A 



. dX 

dgtth _ 
dA 



1 



16tt 2 
1 



2AX 2 -(A-12gl h )X + B 



ugtth 



9tth 

16tt 2 I 12 



\lg' 2 - 27g' 2 - 96g 2 } (123) 



where, 



A = 3(3.g + 9 ) 



B=~-[2g* + tf+g> 2 ) 2 



(124) 



We have omitted the equations for the running of the gauge coupling con- 
stants. They can be found in Rcf. 48. 

In order for electroweak symmetry breaking to occur, the quartic cou- 
pling must remain finite at all scales, 

1 



A(A) 



> 0. 



(125) 



As A is evolved from the weak scale to higher energy, it grows and eventually 
becomes infinite at some large scale. When A becomes large, the scaling of 
Eq. 124 can be approximated, 



A 



dX _ 3 
dK ~ 2^2' 



Eq. 126 is easily solved, 
1 



A(A) 

The position where A (A) 



X{Mh 



3 / A 
-2^1og( — 



(126) 



(127) 



oo is known as the Landau Pole. Since A is 



related to the Higgs mass, M% = 2u 2 A, the requirement that the theory not 
approach the Landau pole gives an upper bound on the Higgs boson mass, 

4ttV 



3 log [A/M h 



(128) 



This means that at the scale A there must be some new physics beyond 
that of the SM which enters. The SM simply makes no sense at energy 
scales larger than A since in this region the Higgs quartic self-coupling is 



c g and g' are the SU(2)l and U(l)y gauge coupling constants. 
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infinite. Assuming no new physics before the GUT scale, A ~ 10 16 GeV, 
gives a bound on the Higgs mass of 49 ' 50 

M h < 160 GeV . (129) 

This bound corresponds to the upper curve in Fig. 26. Non-perturbative 
studies using lattice gauge theory confirm the validity of this limit. 51 

Conversely, we require that the potential be bounded from below, which 
is equivalent to the restriction that A > at large The scaling of the 
quartic coupling near A = is, 

dX 1 



A 



dk 16tt 2 



B - I2g 4 tth 



(130) 



which can be solved, 



A(A) = A(M l) + ^klo g (-A.). (131 , 

The requirement that A(A) > gives a lower bound on the Higgs mass, 
Ml 12gj h -B J A A 

2«2 > l 67r 2 1 ° S {mJ ■ [162) 

A complete approach must include the full set of renormalization group 
equations, including the running of the coupling constants for large values 
of $. 

The restrictions on the Higgs mass from consistency of the theory are 
shown in Fig. 26, where the allowed region for a scale A is the area between 
the curves, known as the "chimney" . At a given scale A, the Higgs mass has 
a lower limit from bounding the potential from below and an upper limit 
from the requirement that Higgs quartic coupling, A, be finite. The exact 
shape of the allowed region depends sensitively on the value of M t . 49 ' 50 

The SM with a Higgs mass in the 130 GeV region is consistent at energy 
scales all the way to the Planck scale. The theory is still unsatisfactory 
however, because it provides no explanation of why /j 2 < 0, as required for 
EWSB. In addition, the Higgs mass is a free parameter in the SM and there 
is no understanding of why it should have any particular value. In the next 
section, we will see that low energy supersymmetry can potentially solve 
both these problems. 



5.2. Supersymmetry and the Top Quark 

TeV scale supersymmetric models explain several mysteries about the SM: 
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Figure26. Limits on the Higgs boson mass as a function of the cutoff scale A. The 
theoretically consistent region lies between the two curves. 49 



• A heavy top quark and mSUGRA (minimal supergravity) explain 
why \i 2 < 0. 

• A heavy top quark and the MSSM (minimal supersymmetric Stan- 
dard Model) explain why the Higgs mass may light, as suggested 
by precision measurements. 



In this section, we give a brief review of phenomenological models of su- 
persymmetry and discuss the relationship between the top quark mass and 
the requirement that supersymmetry exist below the TeV scale. 
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Supcrsymmetry is a symmetry which relates the masses and couplings 
of particles of differing spin. To each SM chiral fermion is associated a com- 
plex scalar and each of the SU(3) x SU (2) l x U(1)y gauge bosons similarly 
acquires a Majorana fermion partner. Supersymmetry connects particles 
of differing spin, but with all other characteristics the same. 52 In an unbro- 
ken supersymmetric model, particles and their superpartners have identical 
masses. It is clear, then, that supersymmetry must be a broken sym- 
metry if it is to be a theory of low energy interactions. There is no scalar 
particle, for example, with the mass and quantum numbers of the electron. 
In fact, there are no candidate supersymmetric scalar partners for any of 
the fermions in the experimentally observed spectrum. The most general 
soft breaking of supersymmetry introduces masses and mixing parameters 
and the theory loses much of its capabilities for specific predictions. 

The Higgs sector in a SUSY model always contains at least two SU(2)l 
doublets. In the supersymmetric extension of the Standard Model, the 
Higgs doublet acquires a SUSY partner which is an SU (2) l doublet of Ma- 
jorana fermion fields, hi (the Higgsinos), which contribute to the triangle 
SU(2)l and U{l)y gauge anomalies as described in Section 2.1. Since the 
fermions of the Standard Model have exactly the right quantum numbers to 
cancel these anomalies, it follows that the contribution from the fermionic 
partner of the Higgs doublet remains uncancelled. 6,7 This contribution must 
be cancelled by adding a second Higgs doublet (and its Majorana fermion 
partner) with precisely the opposite U(l)y quantum numbers from the first 
Higgs doublet. 

All supersymmetric models thus have at least 2-Higgs doublets and 
so the minimal supersymmetric model has 5 physical Higgs bosons: 
h ,H° 7 A° 7 H ± . The ratio of the vacuum expectation values of the neu- 
tral Higgs bosons, 

tan/3= — (133) 

V\ 

is a fundamental parameter of the theory. At tree level, the neutral Higgs 
boson masses are predicted in terms of two parameters, which are usually 
taken to be tan/3 and the mass of the pseudoscalar Higgs boson, Ma- The 
neutral Higgs boson masses are then predicted, 

Ml, H = \{M\ + Ml) T y(Mi + M|)2-4AflM 2 cos2 2/3 (134) 
which gives an upper bound on the lightest Higgs mass: 



M h < M z | cos 2/3 |< M z 



(135) 
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which is excluded experimentally. (In a SUSY model, the LEP and LEPII 
results require M h > 89.8 GeV 4 ) 

The Higgs mass prediction in a supersymmetric model receives large 
radiative corrections from top and stop squark loops. At 1-loop, and ne- 
glecting possible soft SUSY breaking tri-linear scalar mixing, the lightest 
Higgs mass becomes: 

M 2 = Ml+Ml + e 1 | (M 2 + M 2 +£) 2_ 4M 2 M 2 CQS 2 ^ 

, 1/2 

-4e(M| sin 2 + M\ cos 2 0) I (136) 

where, 

3G F M* /M?\ 

e = — =^ ^log —4 , 137 

V2^ 2 sin 2 /3 & \MiJ' K ' 

and M t is the average stop squark mass. For unbroken supersymmetry, the 
top and stop have the same masses, M t = Mt, and the bound of Eq. 135 is 
recovered. The mass splitting between the top and stop in a broken SUSY 
theory changes the upper bound on the Higgs mass to, 

Ml < Mf cos 2 20 + e sin 2 0. (138) 

The large M t 4 corrections raise the Higgs mass bound above the experimen- 
tal limit, as shown in Fig. 27. 

There are many analyses 53 which include a variety of two-loop effects, 
renormalization group effects, etc., but the important point is that for given 
values of tan/3, the scalar mixing parameters, and the squark masses, there 
is an upper bound on the lightest neutral Higgs boson mass. For large 
values of tan the limit is relatively insensitive to the value of tan and 
with a squark mass less than about 1 TeV , the upper limit on the Higgs 
mass is about 110 GeV if mixing in the top squark sector is negligible. For 
large mixing, this limit is raised to around 130 GeV. 



5.3. mSUGRA 

The minimal supersymmetric model (MSSM) is constructed by adding to 
the Lagrangian all of the soft supersymmetry breaking terms allowed by the 
gauge symmetries. These include masses for all the new SUSY particles 
and tri-linear scalar mixing terms. The most general Lagrangian of this 
type has 105 new parameters beyond those of the Standard Model. In 
order to make definite predictions, the simplifying assumption that the soft 
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Figure27. Upper bound on Higgs mass in MSSM at one loop and with no soft trilincar 
mixing. 53 



parameters unify at Mqut is often made. (The resulting model is usually 
called mSUGRA or the CMSSM) . In this framework the model is completely 
specified by five parameters: 



• Common scalar mass at Mqut, "fho 

• Common gaugino mass at Mqut, Mi/ 2 

• 1 Higgs mass-squared, m\ 2 

• 1 tri-linear coupling at Mqut, A \f 

• sign(/i) 
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The parameter fi describes the mixing between the two Higgs doublets. At 
the GUT scale, the Higgs masses are, 

MIh(M gut ) = m 2 + fl 2 , (139) 

while the squark and slepton masses are 

M?(M GUT ) = Mf(M GUT ) = ml . (140) 

The theory must explain why the lightest neutral Higgs mass-squared be- 
comes negative at the weak scale and breaks the electroweak symmetry, 
while the scalar fermion masses remain positive so as to not break color 
and electromagnetism. 

In an mSUGRA model, the masses and parameters are evolved using 
the renormalization group equations (RGE) to scale from Mqut down to 
the weak scale in order to make predictions for the masses at the weak scale. 
The large top quark mass plays an important role here. We can understand 
the basic scaling by considering only the Higgs doublet, H 2l d the scalar 
partners of the third generation SU(2)l fermion doublet, (Jl^l) = Q3, and 
the scalar partner of the right-handed top quark, Ir. The renormalization 
group scaling is given by: 54 

^{Zg 2 tth Xt-Z9 2 M 2 2 -g' 2 Ml) 
-L^Xt-^Ml-fg^M 2 ) 
^( 9 2 th X t - fg 2 M 2 - 3g 2 M 2 \g' 2 M 2 ) 

X t ee M 2 H2 + M 2 3 + M 2 R + A 2 . (141) 

-^3,2,1 are the gaugino masses associated with the SU (3) x SU{2) l x U(1)y 
gauge sector, g U h is the coupling of the Higgs to the top quark in the MSSM, 
and At is the soft SUSY mixing term in the top sector. The QCD term 
(proportional to g 2 s ) makes M 2 R , Mq 3 increase faster than Mjj 2 as the scale 
is reduced from Mqut- The Yukawa term (proportional to g t th) drives 
Mfj 2 < 0, while the other scalar mass-squared terms remain positive. This 
behavior is a consequence of the large M t value. 



dlog(A) 

dM 2 R 
dlog(A) 

dM 2 3 
rflog(A) 



where, 



d H 2 is the SU(2) L doublet with isospin Y = +|. 
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Figurc28. Typical pattern of masses in an mSUGRA model. 55 . 



Neglecting the gauge couplings, we have: 



d 



dlog(A) 



9tth 

8ir 2 



which can be easily solved, 

A4 2 (A) = Mi(MauT) 



M GUT 



(142) 



(143) 



For heavy M t , M% becomes negative and triggers EWSB. It is interesting 
that this behavior requires M t ~ 175 GeV. 55 

The RGE scaling from the GUT scale to the weak scale in mSUGRA 
type models yields a a typical pattern of masses, as seen in Fig. 28. The 
gluino is the heaviest gaugino, while the squarks are significantly heavier 
than the sleptons. 



5.4. Charged Higgs Decays 

The pattern of top quark decays is changed in the MSSM from that of the 
SM, due to the extended Higgs sector. The presence of a charged Higgs 
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boson allows the decay, t — > if the mass of the charged Higgs is less 
than M t —mb- This decay looks similar to t — * Wb and can have a significant 
rate. The coupling of the charged Higgs to the t b system is, 



and so this decay is relevant both for very large and for very small tan /3. 
For small tan/3, H + — > cs and H + t*b, while for large tan/3, 7J + — > tv. 
These decays would suppress the rate for top decays to leptons plus jets in 
the SM top search. In Run I, the Tevatron experiments placed a limit on 
M H+ of roughly M H + > 100 GeV for tan/3 < 1 or tan/3 > 50. 56 

5.5. The Top squark 

The mass-squared matrix for the scalar partners of the left and right handed 
top quarks {Il^r) includes mixing between the partners of the left- and 
right- handed top: 



The mixing is proportional to the top quark mass and so the large value 
of M t can cause one of the top squarks to be relatively light. If the scalar 
mass scale, Mq 3 is much larger than Mz, M t , A t , then the stop masses are 
roughly degenerate ~ Mq 3 . On the other hand if Mq 3 ~ Mz, M t ,A t , then 
the large mixing effects drive the stop mass to a small value and the stop 
squark becomes the lightest squark. 

It is possible to search for a light stop in the decays of the top, t — > tx°, 
where x° is the lightest neutral SUSY particle. The limits then depend on 
the assumed branching ratio for the top into stop decays. The Run I data 
was sensitive to stop and chargino masses in the 100 GeV range. 57 

5.6. The Top Quark and Dynamical Symmetry Breaking 

The top quark has a special role in technicolor models where a it conden- 
sate can play the part of the Higgs boson in generating EWSB. This can 
happen in models where there is a new strong interaction felt by the top 
quark, but not by the lighter quarks. The proto-type model of this type 
is called topcolorJ A review of this class of models can be found in Ref. 



9tH+b ~ m t cot PP + + rut tan /3P_ 



(144) 



/ Mq 3 + M t 2 M t {A t - fi cot 
+ M 2 z L t cos 2/3 

\ M t {A t - fi cot P) + M 2 z R t cos 2/3/ 
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47 . Topcolor models are constructed by expanding the strong gauge group 
from the SU(3) of QCD to, 

SU(3) tc x 5(7(3) -» 5(7(3) QC D (145) 

SU(3)tc couples only to the 3rd generation with a coupling g tc , while the 
second SU(3) gauge symmetry couples only to the first two generations. 
The symmetry is broken to SU(3)qcd at a scale M. Below M, the gauge 
bosons of the SU(3) tc group, the "topgluons", are massive and lead to 
effective 4-Fermi interactions between the top quarks: 

L ~j^(Q3Ll^Q3L)(t R r^-t R ), (146) 

where T A is the SU(3) generator in the adjoint representation and Q%l is 
the left-handed (f^, &l) doublet. If the SU(3)t c coupling, g tc , is larger than 
a critical value, a top condensate forms, (it), which breaks the electroweak 
symmetry and generates a mass for the top quark. One must arrange that 
(it) 7^ and (bb) = 0. This is usually done by adding an additional U(l)h 
gauge symmetry under which the t and b quarks transform differently. 

This class of models can be tested at the Tevatron and the LHC by 
searching for effects of the top gluons and the new Z', and through B 
decays. 59 Since the Z' boson corresponding to the U(l)h and the top gluons 
of SU(3) tc couple differently to the third generation and to the lighter 
quarks, they lead to flavor changing neutral currents. Present limits from 
B meson mixing require that the topgluon be heavier than 3 — 5 TeV and 
M z > > 7 - 10 GeV. 58 

6. Conclusions 

The story of the top quark is far from over. Its properties and its role in 
the SM and in physics beyond the SM will be further elucidated in the 
emerging Run II data at the Tevatron. The increased statistics of Run II 
will enable more precise measurements of the mass, couplings, and decay 
properties of the top. The LHC will provide even more insight into the role 
of the top quark. It is possible that surprises are just around the corner! 
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